Microscopies of disordered two-dimensional electron gases under high magnetic fields: 
Equilibrium properties and dissipation in the hydrodynamic regime 
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We develop in detail a new formalism [as a sequel to the work of T. Champel and S. Florens, 
Phys. Rev. B 75, 245326 (2007)] that is well-suited for treating quantum problems involving 
slowly varying potentials at high magnetic fields in two-dimensional electron gases. For an arbitrary 
smooth potential we show that the electronic Green's function is fully determined by closed recursive 
expressions that take the form of a high magnetic field expansion in powers of the magnetic length 
Ib- For illustration we determine entirely Green's function at order 1%, which is then used to obtain 
quantum expressions for the local charge and current electronic densities at equilibrium. Such results 
are valid at high but finite magnetic fields and for arbitrary temperatures, as they take into account 
Landau level mixing processes and wave function broadening. We also check the accuracy of our 
general functionals against the exact solution of a one-dimensional parabolic confining potential, 
demonstrating the controlled character of the theory to get equilibrium properties. Finally, we show 
that transport in high magnetic fields can be described hydrodynamically by a local equilibrium 
regime and that dissipation mechanisms and quantum tunneling processes are intrinsically included 
at the microscopic level in our high magnetic field theory. We calculate microscopic expressions 
for the local conductivity tensor, which possesses both transverse and longitudinal components, 
providing a microscopic basis for the understanding of dissipative features in quantum Hall systems. 

PACS numbers: 73.43.-f,72.15.Rn,73.50. Jt 



I. INTRODUCTION 



General motivation 



Almost thirty years after the discovery of the quan- 
tum Hall effecfcii 2 -, two-dimensional electron gases under 
magnetic fields continue to attract a considerable inter- 
est both experimentally and theoretically, and have re- 
vealed a rich world of surprising physics. Newly discov- 
ered features concern, e.g., the zero resistance states un- 
der microwave illumination^ and the peculiar Landau- 
level quantization in graphene^. Concerning the integer 
quantum Hall effect (IQHE) itself, direct local imaging 
techniques^ 6 -' 7 -' 8 - have revealed new electron-electron cor- 
relation phenomena and allowed a finer understanding of 
the microscopic ingredients at work. 

On the theoretical side, the quantization of the Hall 
resistance observed in the IQHE relies on the under- 
standing of the quantum motion of charged particles 
in a two-dimensional disordered electrostatic landscape 
in the presence of a strong perpendicular magnetic 
fielc£d°dL±^^iLiL±Li£. As the main effects of the 
electron-electron interaction can be taken into account 
in the integer regime at the single-particle level, using 
a Hartree approach to screenin g 19 : 20 , the calculation of 
equilibrium properties, such as the local electronic den- 
sity and the distribution of permanent currents through- 
out the sample, can be carried out from a one-particle 
random Schrodingcr equation. The precise resolution of 
this problem constitutes the first and important step to- 



ward a microscopic description that underlies the more 
complex nonequilibrium phenomena of the quantum Hall 
effect in its generality. 

Despite the overall good understand- 
ing gained after several decades of 
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pie and general microscopic approach for the physics 
of quantum Hall systems is surprisingly still lacking. 
Computer-based simulations have been propose d 26 ' 27 , 
but rely on heavy numerics in the case of two-dimensional 
disordered potentials, and are not well suited for the 
computation of out-of-equilibrium properties. Even in 
the linear response regime, they are certainly unable to 
address minute aspects such as the tiny deviations to 
the Hall resistance quantization. Analytical approaches 
are better adapted to formulate out-of-equilibrium 
calculations, but face the need to reliably handle in a 
self-consistent screening theory both the formation of 
Landau energy levels at quantizing magnetic fields and 
the complexity of the random potential. 

At present, the theory of the integer quantum Hall ef- 
fect relies on two main cornerstones, which do await a 
unified description. On the one hand, equilibrium den- 
sity profiles arc generally computed within a scmiclassical 
Thomas-Fermi approac h 19 : 20 ' 28 : 29 , leading to a descrip- 
tion of the quantum Hall liquid in terms of compressible 
and incompressible regions. These calculations have how- 
ever proved to be only qualitatively accurate^ -, and fail 
at low temperature, where quantum broadening due to 
the electronic wave function becomes important. Trans- 
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port properties are on the other hand conveniently formu- 
lated in the Buttiker edge state formalis m 19 ' 24 i 28 , which 
nevertheless needs input from more microscopic calcula- 
tions of the bulk properties. This scattering theory be- 
comes also very cumbersome to describe dissipative fea- 
tures of macroscopic samples. An alternative successful 
scmiclassical approach to transpor t 29 ! 30 ' 31 ! 32 ' 33 assumes 
local Ohm's law at a phenomcnological level. A simple 
and controlled quantum approach to both equilibrium 
and out-of-equilibrium properties of quantum Hall fluids 
is thus clearly needed. 



B. Review of the high magnetic field approaches: 
The semiclassical limit 

A popular approach to the IQHE is the high magnetic 
field limit, in which case the center-of-mass motion of 
the electron becomes essentially classical. A quantum 
description is only kept for the orbital effects associ- 
ated with the Landau-level formatio n 10 ' 12 ' 13 ' 14 ' 15 . The 
basic physical idea behind these works is that the ef- 
fective potential seen by the electrons in quantum Hall 
systems is quite smooth at the scale of the magnetic 
length Ib = \Jhc/\e\B (here h is Planck's constant di- 
vided by 2ir, c is the velocity of light, B is the magnetic 
field strength, and |e| is the absolute value of the electric 
charge). This permits a simple mathematical treatment 
of the Schrodingcr equation using as a small parame- 
ter the ratio of the magnetic length to the typical cor- 
relation length of the random potential. This point of 
view is certainly vindicated experimentally by the fact 
that the impurities in semiconducting heterostructures 
are located outside the two-dimensional layer of conduc- 
tion electrons, while Ib is an extremely small length scale 
which falls below 10 nm for magnetic fields above 5 T. 
It is tempting to believe that all aspects of the quantum 
Hall effect should be captured accurately in this limit. 

While the idea is certainly not new, it is interesting to 
note that no fully quantum treatment of the high field 
regime is currently available, which is the issue we want 
to address in this paper. Focusing first on the scmiclas- 
sical corrections to the B = oo limit, it is known that 
systematic calculations are quite cumbersome, even at 
the lowest orders in the Ib expansion, due to the Landau 
level mixin g 13 ' 34 , so that new tools are certainly needed. 
A first technical step in this direction was made by two 
of us in a recent publication^ 5 , by introducing well-suited 
coherent states Green's functions. These so-called vortex 
states with the quantum numbers v = (m, R), where m is 
the Landau-level index, and R the position of a localised 
vortexlike wave function, form an overcomplete basis of 
eigenfunctions with no preferred symmetry, in contrast 
to the widely used translation-invariant Landau states or 
the rotation-invariant circular states. They thus permit 
a great adaptability to the spatial variations of the lo- 
cal electric fields, coming from cither random impurity 
donors, confinement potentials due to external gates, or 



macroscopic voltage drops. This formulation also allows 
one to classify Landau level mixing processes in a simple 
and natural manner, an important point for the inves- 
tigation of quantum transport properties, as the matrix 
elements of the current density necessarily relate adjacent 
Landau levels. 

Our first implementation of this technique^ 5 - has 
demonstrated, not surprisingly, that the usual semiclassi- 
cal approach to the quantum Hall effect (limited to spec- 
tral properties) could be easily recovered by a straightfor- 
ward expansion of the vortex Green's function in powers 
of the magnetic length. In this view, the vortex state 
coordinate R can be identified in the limit Ib — > with 
the slow classical center-of-mass drifting motion, which 
completely decouples from the faster cyclotron motion. 



C. Toward a unified quantum description at high 
magnetic field 

The present paper has two aims. First, we want to 
provide an accurate quantum treatment of the local equi- 
librium properties of quantum Hall systems. For this 
purpose, we offer simple functionals of the arbitrary lo- 
cal electrostatic potential that describe both the local 
charge and current densities. These results may have 
important bearings for microscopic modelizations of real 
devices based on Hartree-Fock or more refined local den- 
sity approximation (LDA) calculations, as they avoid the 
numerical costs in solving the random Schrodingcr equa- 
tion in a magnetic field. The knowledge of the current 
density functional can be used in a second step to ob- 
tain out of equilibrium transport equations, which take a 
simple hydrodynamic form at high magnetic fields. This 
step allows us to derive microscopically a simple and lo- 
cal expression for the conductivity tensor. We show that 
in contrast to the well-known drift contribution to the 
transverse Hall conductivity^, dissipative longitudinal 
components first appear at order l B . These contributions 
had not been obtained previously in the literature to our 
knowledge. This method also allows one in principle to 
derive microscopically the dominant nonlocal corrections 
to local Ohm's law. A general understanding of the dissi- 
pative features in the integer regime seems now possible 
at the microscopic level. 



D. Organization of the paper 

Because the paper involves several novel technical de- 
velopments, we hereafter guide the reader through the 
main results obtained. Sec. [TT]is used to introduce vortex 
Green's functions and reformulate in a more systematic 
manner the results obtained in Ref. l35l . Two important 
formulas are found that determine completely Green's 
functions both in the vortex coordinates [Eq. (p?o| ] and 
in the electronic coordinates [Eq. (|44p]. These are the 
starting point for the computation of all physical obscrv- 
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ables. The local electronic charge density is thus derived 
up to order l 2 B in Sec. IIIII and is given by formulas (f5U)) . 
(|56|) and (|58|) . We emphasize beforehand that all these 
expressions take into account quantum smearing effects 
from the wave function and extend the semiclassical re- 
sults (also derived in this section) to much lower tem- 
peratures. Similarly, the (equilibrium) local electronic 
current density is computed to the same order in Sec. 
HVl and is given by Eqs. ([72]), ([C6]), ([C8]), and 

(|C10p . Again, the semiclassical current density can be 
obtained from these expressions and is given by Eqs. (|73|) . 
(|75|. and (fCTT]) - (fCT4|) . We then provide in Sec. |V]two 
important checks of our theory against an exactly solv- 
able model of a one-dimensional parabolic confinement 
potential. First, the analytic semiclassical formulas for 
the local observables obtained in Sees. IIIII and IIVI are 
compared with the strict expansion in / b from the exact 
model and shown to match precisely, strengthening the 
mathematical foundation of our theory Second, a more 
quantitative comparison is made between the quantum 
expressions obtained for the electronic density in Sec. 
IIIII and the exact results at finite values of the magnetic 
length. This shows that the expansion proposed here 
is converging quite rapidly, even for a confining poten- 
tial that is not exceedingly smooth. Because the vor- 
tex states do not favor any special symmetric situation, 
similar quantitative results should be obtained for an ar- 
bitrary two-dimensional smooth disordered potential. In 
the limit of zero temperature, this comparison also shows 
the need for a resummation of the quantum expressions, 
to infinite order in Ib, modifying both the vortex wave 
functions and energies, and allowing a possible connec- 
tion to the edge state picture. Finally, Sec. |VT] investi- 
gates nonequilibrium properties in the integer quantum 
Hall regime and provides a microscopic derivation of the 
local conductivity tensor in the semiclassical regime [for- 
mula (| 108[) ] . The origin of dissipation is discussed, and a 
general conclusion showing future directions of our work 
closes the paper. Some extra technical details are given 
in several appendixes. 



z direction). In the symmetrical gauge 



he l 



A ( r ) = o^ 



e| 2Z| \ 



the vortex wave functions, with quantum number v = 
(to,R), are written in terms of the electronic variables r 
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*m,R(r) = 

x exp 



1 



x - X + i(y - Y) 

\/2Trm\lB V V%Ib 
(x - X) 2 + {y- Y f + 2i(yX - xY) 



4/2 



The associated energy levels read 

Em,n = (m+ 1/2) hu c = E n 



(2) 



(3) 



where lu c = \e\B/m*c is the cylclotron pulsation. The 
energy levels are independent of the position R, and 
their quantization is uniquely related to the (topological) 
quantization of the circulation of any paths enclosing the 
position r = R which corresponds to a phase singularity 
of the wave function \I/ m .R,(r) (note that this wave func- 
tion vanishes only at the point r = R where its phase is 
ill-defined: It describes a vortex). 

As is clear from Eq. the Landau energy levels are 
highly degenerate, so that there is a great freedom in the 
choice of a basis of states. However, a judicious choice 
for the set of quantum numbers appears essential when 
considering perturbations that lift this huge energy de- 
generacy. A peculiarity of the vortex states, which could 
appear at the first glance as a drawback, is that they 
are nonorthogonal with respect to the degeneracy quan- 
tum number R. Indeed, the overlap between two vortex 
states is given by 

(viWi) = <W,m 2 (R1IR2) (4) 
(Ri - R 2 ) 2 - 2ii ■ (Ri x R 2 ) 



(R1IR2) = exp 



41% 



(5) 



II. HIGH FIELD EXPANSION WITHIN THE 
VORTEX STATES REPRESENTATION 

A. Vortex states 

The vortex states under which our quantum high mag- 
netic field theory reposes are eigenstates of the free 
Hamiltonian 



H 



1 



2m* 



e \ 2 
-A(r) 

c 



(1) 



describing a single electron of effective mass m* and of 
charge e = — |e| confined in a (xy) two-dimensional plane 
under a perpendicular magnetic field B (pointing in the 



where z is the unit vector along the perpendicular mag- 
netic field. 

On the contrary, the other well-known eigenstates of 
Hq, the Landau and circular basis states, which are com- 
monly used for quantum calculations, are orthogonal. 
But since they are highly symmetric states, they lead to 
unsolved technical difficulties when considering a random 
potential in high magnetic fields, which mixes in a very 
complicated way the two quantum numbers. The vortex 
basis, which has no intrinsic symmetry (the nonorthog- 
onality of the vortex states arises from this property), 
allows one to overcome this drawback. The possibility^ 
to work with this basis is, in fact, provided by the co- 
herent character of the vortex position degree of freedom 
(the algebra obeyed by the degeneracy quantum number 
R is that of coherent states). 
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B. Dyson equation in the vortex representation 

From now on and throughout the paper, we consider 
that the Hamiltonian contains in addition to the kinetic 
part Hq a potential energy term V(r), which we let com- 
pletely unspecified 



H = H + V(x,y). 



(6) 



The Dyson equation written within the vortex represen- 
tation \v) = \m, R) then takes the forrn^ 

(w - E mi ± i5)G%f U2 (u;) = {v x \y 2 ) + £ V w G^(u), 

"3 

(7) 

where G^^ 2 (ui) are retarded and advanced Green's func- 
tions connecting two vortex states v\ and vi (in the en- 
ergy representation). The sum over the vortex quantum 
numbers v appearing into the Dyson equation stands for 



E = E 



dR 

m=0' J B 



(8) 



The matrix elements of the potential V(r) in the vortex 
basis are given by 



J d 2 r^(r)*; ii Ri (r) * m2jR2 (r) 



(Ri|R 2 )^ i;i 



(9) 



where, for a practical purpose which will appear obvi- 
ous in the following, the overlap between the two vortex 
states has been extracted. Similarly for retarded and ad- 
vanced Green's function, we extract the vortices overlap 



(10) 



where the dependence on frequency u is not explicited 
anymore, in order not to burden the expressions. Sub- 
stituting expressions © and (fTU)) in Eq. Q, we get a 
Dyson equation for the function g^'^ 2 which reads 



(w - E mi ± i6)g%£ 2 

where 
(R 1 |R3)(R 3 |R 2 ) 



a R ' A 



^ (R 1 |R 3 )(R 3 |R 2 ) 



<Ri|R 2 



exp 



(R 3 



(c 12 + id u x z)Y 



2/ 2 



(12) 

We have introduced here the center-of-mass coordinates 
C12 = (Ri + R 2 )/2 and the relative coordinates di 2 = 
(R 2 - Rx)/2. 

Provided that V(x, y) is an analytic function of both x 
and y, the reduced matrix element v Vl - V2 of the potential 



appearing in Eq. ([9]) can be written^ as a series in powers 
of the magnetic length I b 



E 

3 

E 



; \ ■? 



s/2J VxW 
(mi + fc)! 5 mi+ k, m2 +j-k 



(13) 



Vl ' va ^ V^W k\{j-k)\ 



x(d x + l d Y ) k {d x - i d Y y- k ^(R)| ci2+ , di2X2 . 

(14) 

Solving exactly Dyson equation [Eq. (fTTj) ] for an ar- 
bitrary potential V is certainly a formidable task. Re- 
markably, however, from the structure of this equation, 
one can show that the function g Ul -u 2 depends on the two 
vortex coordinates Ri and R 2 through the special com- 
bination ci2 + idi2 x z only. Indeed, let us differentiate 
Eq. flT]) with respect to the first vortex position, 



{u-E mi ±iS)(d Xl -id Yl )g^ 2 

E/o . a \f R A (Rl|R3)(R3|R2) 

/3 " {^9^ 2 (Ri|R2) 

Then, noting that from Eqs. ((T2J) and fH]) 



(15) 



. (Ri|R 3 )(R 3 |R 2/ 
(d Xl - id Yl ) < , p r > = 0, (16) 



(Ri|R 2 ) 

(d Xl - id Yl )v„ llU3 = 0, (17) 



and considering Eq. (|15p . we arrive to the relation 

dx 1 g%£ a =idy 1 g%£ a . (18) 

We can establish similarly from the other Dyson equation 
(i.e., G = G Q + GVG ) that 



dx 2 g^- U2 = -idy 2 g p 



(19) 



We thus deduce from these two relations [Eqs. (Tj"8")) and 
(p~9|) ] that the function g depends on the vortex positions 
Ri and R 2 in the following way: 

^ R 1 +R 2 + t (R 2 -R 1 )xz ^ 

9vi\v 2 — gm 1 ;m 2 I _ I • (.^"j 



This exact result implies that vortex Green's functions 
will be entirely determined once the function g Vl -,v 2 at 
coinciding vortex positions Ri = R 2 = R are known 
(provided it is analytic in the complex plane). This task 
is addressed in Sec. II C. 



C. High magnetic field expansion of vortex Green's 
function 



We are mainly interested in the high mag netic fie ld 
regime, i.e., when the magnetic length lg — yJhc/\e\B is 
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small compared to the typical length scale of the (pos- 
sibly random) potential V(r). We aim at solving the 
Dyson equation [Eq. () 1 1 [) ] as a systematic expansion in 
powers of Ib, i-e., expanding the function g Vl -^ 2 as 

+°° / 7 \ j 
=£(-%) (21) 

This expansion is possible because, using the change in 
function (|10[) . the nonanalytic dependence on the mag- 
netic length Ib which was contained in the first term of 
the right-hand side (rhs) of Eq. ([7]) has been fully trans- 
ferred to the overlap terms [Eq. (p~2|) ] appearing in the 



integral contribution of the Dyson equation [Eq. 
At large magnetic fields (i.e., small Ib) and when Ri is 
close to R.2 the main contribution to the integral over 
R3 in Eq. (|lip comes when R3 is near both positions 
Ri and R2. Because Green's function (|20|) depends on 
a linear combination of the two vortex locations, it is 
enough to calculate vortex Green's function at coincid- 
ing points Ri = R2 = R, so that, from Eq. (fTTj) . only 
the value of the function u mi) B,- m8 ,Rj,ff m3 ,Rj; m j I R has to 
be considered. The Dyson equation can now be solved 
by expanding the nonlocal functions g and v around co- 
inciding points using a Taylor series at R3 close to R, 



<?>T13,R.3; 



^ [(X 3 -X)-i(Y 3 -Y)] k k 



"mi ,R.;m3 ,R,3 



k=0 
+00 

E 

fc'=0 



[(X 3 -X) + i(Y 3 -Y)] k v m, 
(d x - idy) u miiTn ,(R), 



(22) 
(23) 



where we have taken into account the spatial dependences of v [see Eq. (fl4|)] and g [see Eq. (|20l) ]. The integral over 
the vortex position R3 in Eq. can then be evaluated using the following property of Gaussian integrals: 



2tt/ 



2 ^mi,R;m3,R3^m3,R 3 ;m2,R e 



(R, 3 -R)^ 



2ttH 



[(x 3 -x)- z(y 3 - Y)] k {(x 3 -x) + z(y 3 - Y)f 



k.k'=0 



(R 3 -R) 



k\k'\2 k + k 



- (d x - idy) fmi;m3 

(R) {d x + idy) 

#m 3 ;m 2 

(R) (24) 



+OO 



k=0 



2 k 



Combining the different series expansions of the matrix elements of the potential [Eqs. (|13|) - (p^|) ] . of vortex Green's 
function g R ' A [Eq. I|2ip]. and of the integral term [Eq. I|25p] in the Dyson equation, the functions gf%£ 2 at coinciding 
points Ri = R2 arc then entirely determined order by order in powers of the magnetic length Ib- In fact, vortex 



Green's function at order l 7 ^ is related to the terms g^ (with I < n) through a closed-ft 



orm recursive 



relation: 



n-l n~l (n-Q/2 



mx+j 



iU( R ) = jSUMEE E £ (^+z9 y )%« ;m2 (R)(ax- l 9y) fc 4i ) i;m3 (R) I (26) 



1=0 j=0 k=0 



m 3 =m 1 -j 
I 



points 



where the function is given by Eq. (fT4|) at coinciding magnetic length. Even for a disordered potential that is 

smooth on the scale of Ib, the question of the accuracy 
and convergence of this expansion has to be addressed. 
We refer the reader both to a general discussion of this 
important point in Sec. Ill Fl and to a concrete comparison 
with an exactly solvable model in Sec. [V] 



v {i) (R) = N 



k=0 



(mi + k)\ 

Vmi!m 2 ! kl(j — k)\ 



x(d x +idy)*(d x ~ idYy-«V(R), (27) 

and the zeroth order contribution g^\ that suffices to 
determine the whole series, is given in Eq. (|29[) below. 

Obviously, the present method generates a systematic 
expansion for vortex Green's functions in series of the 
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D. Vortex Green's functions up to order 1% 

For the calculations to follow in the rest of the paper, 
vortex Green's functions up to order lg will be needed, 
and these useful expressions are given here. At leading 
order (zeroth order in magnetic length) the equation de- 
termining the function is trivially found by setting 
k = in formula (|25j) . which is then reported in Dyson 
equation [Eq. fTTjl]. 

(u-E mi ±i%W. ro2 (R) = S mum2 + V(R)g£\. m2 (R). 

(28) 

This equation is entirely closed and yields straightfor- 
wardly 

with £,„(R) = E m + V(R). Green's function at leading 
order is diagonal with respect to the vortex circulation 
quantum number m. We regard this robustness of m 
independently of the detailed form and strength of the 
potential V as a signature of its topological nature. We 
see that, in addition to a kinetic term (E m ), the energy of 
the vortex state £ m (R) now also contains the value of the 
potential energy V(R) at the vortex location, which lifts 



the huge degeneracy of the Landau levels. This leading 
order of the calculation clearly corresponds to the strict 
semiclassical limiti^iiiiH^i at lg = 0. 

All subleading contributions are straightforwardly de- 
termined using the recursive relation (|26[) . which for 
n = l gives the order Is contribution 

S^( R ) = 5^; mi ( R ).9^; m2 ( R )^;™ 2 ( R ) (3°) 

^mi ;rri2 (^) 

; (w - (R) ± iS) (w - U 2 (BO ± iS) ' 

(31) 

where from equation (1141) 

{Ox - idy)] V(R). 

(32) 

We thus see that a mixing between adjacent Landau lev- 
els appears in the presence of a gradient of the potential 
V. 

For the determination of the function g^ 2 \ the matrix 
elements of the potential at order l 2 B are needed, which 
read from (fT4|) 



v {2) (R) = - 



1 max (mi, 7712)! 



2 \Jm\\mi\ 
Recursion relation (|26[) at this order gives 



(dx+idy) + (d x - idy) + 2(mi + l)<5 mi , m2 A R V(R). (33) 



o (2) (R) = o (0) 



^ ; ™ 2 (R)^; m2 ( R ) + E W ™( R )^;™ 2 ( R ) + <5 mi ,rn 2 V R y(R) • V R5 W. mi (R) 



(34) 



The function g( 2 ) consequently contains diagonal elements (mi = 7712) and elements mixing Landau levels separated 
by an energy of 2hu> c (terms with mi = m2 ± 2), 



7 (2) (R) - 8 

Jm 1 ;m 2 \ CK ) ~ "nil, 



(mi + 1)— 2— + + + j 

(a x + *9y) 2 ^ , Px + *9y)^] 5 



+Vmi +Tvmi+2 <5 mi+2 ,. 



+Vm 2 + lVm 2 + 2<5 mi , m2+2 



(5 A -- J 9y) 2 F , Px-i9y)l/] 2 



(35) 



where we have introduced the short-hand notation Lu m = ncnts that are needed for subsequent calculations. 
lu - £ TO (R) ± id. 

As for third-order Green's function, we shall not write 
here the full expression, which is rather cumbersome. 
The derivation of these terms from Eq. (|26p is however 
straightforward, and Appendix [X] provides the compo- 
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E. Green's functions in the electronic 
representation 



The aim of this section is to connect local vortex 
Green's function, determined previously in the magnetic 
length expansion, to physical observablcs. For this pur- 
pose, we need to express Green's functions in terms of the 
electronic positions r, which, thanks to the completeness 
relation satisfied by the vortex states 3 ^, can be obtained 
as 

G(r,r» = £ G VUV2 (w)** 2 (r')*^ (r). (36) 



Rewriting in terms of the vortex location R and circula- 
tion m, this expression reads 



H,Ri;m 2 ,R 2 — 3mi;m 2 ( c 12 + *d 12 X z) 

= g r mi;m2 (R + 2idi2 x z) 



(40) 
(41) 



\- [2i(di 2 x z) ■ V R ] 3 

2^ J, Sm i; m 2 (R) (42) 



J' 



3=0 

y> y-y (di2x - idi2y) k (di2x + idi2 V ) 3 ~ k 

j=0 k=0 



k\(j - k)\(-l) k 

x(dx + id Y ) k {dx - td Y y- k g mi , m2 (R) 

(43) 



where we have used 2i(di 2 xz)-VR = (di2 X +idi2 y )(dx + 
idY) — (di2 X — idi2y){dx—idY) and then applied the bino- 
mial theorem. Inserting expansion (|43]) into Eq. ([37]) and 
using Eqs. |38|) -([39 |) . we can then perform the integral 
over the relative coordinates di 2 to finally obtain 



G(r,r',w) — f— „ f— „ G mi R i;m2 r 2 (cj) , Ai 2 R ^ , 

J 2nl B J 2tt1 2 b G(r,r ,u) = ^ tf m /, R (r )* TO , R (r) 

B rn to ' 



(37) 



+ DC 



(~ A, 



Besides the double integral in the above formula, the 
difficulty we immediately encounter is that nonlocal vor- 
tex Green's function is in principle needed. Again, we are 
going to see that the key formula (f2"0"|) allows us to refor- 
mulate this expression in terms of local vortex Green's 
function determined in Sec. Ill Gl Inserting expression 
(|10[) . we first write the different exponential factors ap- 
pearing into the integrand of the expression ([3T]) as 



k=0 



F. On the convergence of the Ib expansion 



(r-Eif-2i(ixRi) : 



(r'-R 2 ) 2 + 2 i(r' X R 2 )-: 



(Ri|R 2 )e 



e 

- 2 i(rxR)-: 



b e 



(r'-R) 2 + 2 i(r'xR) 

7Ti 



,(38) 



where R = C12 — «di2 x z is a complex combination of 
the center-of-mass and of the relative vortex coordinates. 
Similarly, the polynomial parts of the vortex wave func- 
tions can be written as 

[x> -X 2 - i(y' - Y 2 )] m2 [x-X 1+ i(y - Fx)]" 11 = 

[x'-X- i(y' - Y)] m2 [ x -X + i(y- Y)] mi ■ (39) 

It thus seems natural to introduce the change in vari- 
ables (Ri,R2) — > (R, di2). The variables X and Y lie 
a priori on lines in the complex plane as a result of the 
complex shift (— idi2 x z). Using the analycity property 
of the functions in the integrand, the contours of integra- 
tion can be deformed to the real axes. The dependences 
on the variables R and d 12 in the function 3mi,Ri;m 2 ,R 2 
arc made separable using Eq. (|20[) and expanding the 
nonlocal g function as 



Equation (|4"4"|) above is clearly remarkable as it con- 
nects local vortex Green's function <? m;m '(R) to the non- 
local electronic propagator G(r, r',w), from which all 
equilibrium physical properties can be obtained. Because 
the vortex wave functions appearing in this expression 
have a finite extension in space of order Ib, the combi- 
nation of Eq. (|4"4"|) with recursion relation (f2"6"|) , which en- 
codes the small Ib expansion of g m - m i(R), allows one to 
systematically obtain quantum expressions for the physi- 
cal observables, i.e., that are naively valid at small but fi- 
nite magnetic length. In contrast, the usual semiclassical 
cxpansio n 13 i 34 i 36 is formulated in a strict Ib — > limit, 
which would appear in our formalism as a further expan- 
sion in powers of Ib of the wave functions in Eq. (|44|) . 
The latter semiclassical expansion, which is analyzed in 
detail in Sees IIII Dl and IIVD1 is clearly asymptotic in 
nature and certainly fails to be accurate at low temper- 
ature, where quantum effects set in (this is explicitely 
demonstrated in Sec. [V] with the comparison to an ex- 
actly solvable model). 

A central question is whether our expansion, per- 
formed order by order in powers of I b for vortex Green's 
function <?(R), does fully capture the quantum effects 
that survive at small but nonzero magnetic length. As 
discussed by several author a 10 ' 12 , the Schrodinger equa- 
tion becomes integrable in this limit, with constants of 
motion associated to equipotcntial lines of the random 
potential. The wave functions are then localized on a 
scale Ib transverse to these equipotcntials and spread on 
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the whole constant energy contour. While this picture is 
certainly appealing, difficulties arise for generating a sys- 
tematic expansion that takes into account Landau level 
mixing. In contrast, the vortex states are flexible enough 
to capture these important contributions but correspond 
to a starting point where wave functions are fully local- 
ized at the scale Ib- Indeed, the first term in our expan- 
sion [Eq. (f2"9")) ] indicates that eigenstates in the B = oo 
limit correspond to equipotcntial points, not lines, and 
this is related to the fact that the vortex wave functions, 
which form an ovcrcomplete basis, become pointlike and 
thus orthogonal in this limit. Physically, one expects that 
quantum fluctuations will play a crucial role as soon as B 
is finite by selecting orthogonal, and therefore more ex- 
tended, wave packets. Mathematically, this phenomenon 
is reflected in our formalism by the presence of terms that 
have to be kept at each order of the small I b expansion for 
vortex Green's functions [e.g., at order l B , these are the 
terms with a single frequency to rn in Eq. (|35p ]. The need 
for resumming this expansion comes in fact from the Tay- 
lor expansion to finite order of nonlocal vortex Green's 
function [Eqs. (|22|) and f|43J)] . Fortunately, as terms at 
an arbitrary order in Ib can be generated through re- 
lation (|26[) . it is possible to achieve a resummation of 
the leading contributions to vortex Green's function. Al- 
though these considerations are beyond the scope of the 
present paper, a first step in this direction is presented 
for the case of the local electronic density in Sec. IV Dl 



III. ELECTRON DENSITY 

A. General expression 

Vortex Green's functions being determined, one can 
then derive quantum microscopic expressions for the lo- 
cal physical observables. The equilibrium local density is 
related to distribution (lesser component) Green's func- 
tion G < in the electronic representation by the general 
formula 



P(r) = ~i J^G<(r,r,u) 
[duj [d 2 H 



(45) 



+°° i / ;2 \ k 

x E.T -f Ar ^'( R )' < 46 ) 



where Eq. (|44j) has been used. The distribution function 
in the vortex basis reads at equilibrium 

-fcmjm'W = m FH« m ,(R)-ff^ ;m ,(R))(47) 

f 



1 + exp [(w - n*)/T\ ' 



(48) 



where the Fcrmi-Dirac distribution function has been in- 
troduced, with jU*/e = <f>o the electrochemical potential 



(which is constant in space at thermodynamic equilib- 
rium). 

The computation of the electronic density at an ar- 
bitrary order in the magnetic length expansion is now 
straightforward using Eq. (|26|) to generate successive con- 
tributions to <?™. m '(R)- An important remark is however 
in order. While equation (|4"6")l involves local Green's func- 
tion relative to the vortex position R, it takes into ac- 
count all Landau level mixing processes (terms with m ^ 
to'). As we will discover in the following calculations, the 
combination of vortex wave functions ty* m , R (r)>I' mi R.(r) 

involves an extra power l B m ' . For this reason, the 
contribution from which couples adjacent Landau 
levels, is actually of order l 2 B and not Ib- Similarly, the 
diagonal (to = m') terms in are indeed of order l B , 
while contributions with to = m! ± 2 [see Eq. (|35|)] are 
overally of order l B and will be discarded in the following. 



B. Electron density at leading order 

At leading order of the expansion in Ib , vortex Green's 
function is given by Eq. (|29p so that the distribution func- 
tion reads 

" ^m?m'( R ) = 27m F (o;)cW S(u - £ m (R)). (49) 

Inserting this in Eq. (|46[) and performing the frequency 
sum, we obtain the local electron density 



p(°)(r) 
l* m ,R(r)| 2 



j2t> +°° 

^T^^mMl'nF (MR)) (50) 



1 



2nm\l% 



R 



V2h 



2m 

exp 



(R 



21% 



(51) 



We note that this zeroth order contribution (|50|) is al- 
ready more powerful than the expression for the electron 
density that is obtained in the strict limit Ib — > of 
infinite magnetic field, 



p(°)(r) 



1 



+00 



2nl 2 R ^ 



nF (£m(r)) ■ 



(52) 



This semiclassical result (f52|) has been widely used in 
the literatur o 19 i 36 as a basis to screening calculations. It 
however ignores the fact that the physical density cannot 
vary faster than the scale Ib as is clear from Eqs. ([50]) 
and (|51| . and this leads to important quantitative dif- 
ferences, especially at low temperatures where quantum 
smearing effects supersede the thermal broadening of the 
density^. 

Thus, expression (|50[) clearly includes important re- 
summations of a purely semiclassical expansion of the 
physical density such as Eq. (|52p . which is naturally en- 
coded order by order in the expansion of vortex Green's 
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functions. Before addressing the question of the conver- 
gence of both types of calculations in Secs lIII Dl and lTVDl 
we compute now the next order contribution to the den- 
sity. 



C. Electron density at order 1% 



As mentioned above, the contribution of order l 2 B to the 
density comes from three origins: The non-diagonal part 



J 



of g( x \ the diagonal part of g^ 2 \ and the term Ar^ ) in 
expression (|46j) appearing with k = 1 . Let us investigate 
these different contributions in turn. 



1. Contribution from g^ 

The contribution from g^ is obtained by inserting 
Eq. (f5T|) in Eq. (j!?Tj) , and reporting vortex Green's func- 
tion in Eq. (f4H|) : 



pW(r) 



M2t> +°° ; 

J^jT £ -jf 1 ™^ 1 [*WW*^W(^ y + WyV ^> + *m,R(r)*m+i,R(r)(9 X V - id Y V)] 



B m=0 

[»F (£m+l( R )) - (£m(R))] 



(53) 



It is useful to note the relation, proved in the Appendix (|53[) can be written in the equivalent form 

m 

(% \ _ _ ^?MUi(R))-nF(( m (R))] 
P W — /o_;2 / / 



Vm+14; +1|R (r)t m , R (r) - 

, m 
>>B / n . a x NT. / \|2 



2tt/ 2 r - 

« m=0 



~(^-»^)x;i*p,R( r )r 

* p=0 

which shows that the product of wave functions 
i i 
c 

sity is indeed of order l B . Using Eq. ([5"3|). the expression we end up with 



1 2 b VrV ■ £ |* Pl R(r) 

p=0 



(55) 



with adjacent Landau indices generates terms that be- Performing one of the discrete sums and an integration 
have as Ib, so that the contribution from g^> to the den- parts and noting that V r |^ m ,R(r)| = - Vr | v l'm,R(r)| , 



P (1) (r) = 



j^T £l*m,R(r)| 2 £ Kf(^(R)) 



- n F (4m (R)) ■ 



huj r 



fkU r 



(56) 



Contribution from 

From Eq. (|46|) the remaining contributions to the electron density at order l B are clearly 



,(2) 



(r) = - J /S/SE^K(r)* m ,R(r) 



r |^ 2 ;l1(R)-|A R ^)<(R) 



Using results (f!29"j) and (f3"5| . we obtain 

72 T3 +°° 



(2) 



^ElMr)| 2 f <(^ m (R)) 



s m=0 
2 



mA R U 



iv R ur 



n'lr (€m(R)) 



iv R ur 



K m + X ) nF (£™+i( R )) - ( 2m + !) n ^ (MR)) + ranjr (£ m _l(R))] 



(57) 



(58) 



The components <5 m ,m'±2 of the function g^ [see Eq. they generate corrections to the density of the order I 
(|35p] have not been included in the calculations since 



B ■ 



10 



The final results for the electronic density up to order 
l 2 B , given in formulas (f50")l . ([55)1 and ([55]) , will be exploited 
in detail by a comparison with an exactly solvable model 
in Sec. [V] In anticipation to Sec. IV Dl we note that all 
these expressions require at very low temperature a re- 
summation procedure, which leads to define renormalized 
energies and wave functions. 

D. Semiclassical density: The strict Ib — > 
expansion 

As already mentioned above, it is possible to express 
the electron density under the form of a strict expansion 



in powers of the magnetic length. This corresponds ex- 
actly to a systematic semiclassical expansion with respect 
to the center-of-mass motion (the orbital motion giving 
rise to the Landau levels is always treated quantum me- 
chanically) . 

In this section, we write down explicitly the first cor- 
rections to the well-known semiclassical expression, i.e., 
Eq. (f52")l . for the electron density. Nonlocal expressions 
(f50")l . ([55]) and ([55]) can be transformed into local ones, in 
a similar way as has been done at the level of the Dyson 
equation [Eq. (fTTj) ]. by replacing any integral over the 
vortex position R in the following way: 



n n n ~\~ OO 

/d 2 R|* m , R (r)| 2 /(R) = /rf 2 R|* m , R (0)| 2 /(r + R) = /d 2 R|* m ,R(0)| 2 J]. 

J J J j = Q 



(R • V r ) J 



»( m + i )!/Z| w 



E 

3=0 



(j!) 2 m! V 2 



^A r /(r), 



(59) 
(60) 



where /(R) represents an arbitrary function of the vortex position. Using Eq. (foT)]) in Eq. (f5U|) . we get the semiclassical 
contribution of order l 2 B arising from the expansion of p^°\ 



B m=0 



n F (f m (r)) + y (m + l)A r n F (f m (r)) 



(61) 



At order l B , the contribution from [Eq. (|5"6"|) ] is readily obtained as 



p«(r) = 



1 + °° / 



(6n(r)) 



+ n F (£ m (r))A r V 



(62) 



Similarly, the contribution from [Eq. I|58p] reads after simplification 



hoo , 2 



2tt/ 2 ^ 2 



hlOr 



- 4 (^m(r)) 



|v r yr 



(63) 



Collecting all the different contributions, we find that the total electronic density in high magnetic fields including 
the first quantum corrections of order l B is given by 



m=0 



Ttt = 



-n F (£ m (r)) -jj— - + j 



»F (fm(r)) + | 



4 (^ m (r))|V r Vr 



In order to physically interpret this result, we can alternatively write 



* +oo 

P^) = ^jT E 



A r U 



nf(fm(r)) + ?|n F (£ m (r)) -r±— + ^ (£ m (r)) 



V r V\ 2 I 2 



-£lm + -)A r n F (U(r)) 



(64) 



, (65) 



where corrections proportional to n' F in p' 2 ' have been absorbed into renormalized Landau-level energies as 



£m(r)=£m(r) + %(m+i) A r F- 



2 ftu; r 



(66) 
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These semiclassical energies [Eq. ([66)) ] have been previ- 
ously found in the literatur o 13 i 34 with techniques based 
on effective Hamiltonians, which neglect Landau level 
mixing, and thus do not allow one to compute the full 
local-density expression (|65[) . The second term appear- 
ing in the rhs of Eq. (f65|) . proportional to l B nF^rV/hoj c , 
reflects the small but nonzero compressibility of the elec- 
tron gas even in the absence of electron-electron inter- 
action. This term and also the third term in Eq. (|65|) . 
both derived from stem from adjacent Landau level 
mixing processes and can be interpreted as small correc- 
tions to the wave function for a smooth potential V(r). 
Clearly the fourth term in Eq. (f65|) . proportional to 
l B A r nF, is only a small correction if the electronic den- 
sity is also smooth at the scale of Ib- This indicates that 
the semiclassical picture breaks down at low tempera- 
ture. In this case, one has to resort to fully quantum 
expressions such as Eqs. (f5U|) , ([56]) and ([58"]) . as will be 
discussed in Sec. IV CI 



calculation o 19 i 20 i 37 . The physical picture that emerged 
from these studies is that the sample separates into either 
compressible regions, where screening of the bare poten- 
tial is almost perfect and the electronic density varies spa- 
tially, or into incompressible regions, where the density 
is almost exactly pinned and the gradient of the effective 
potential is nonzero. Further work 3 ^ has however shown 
that important deviations result from a better resolu- 
tion of the self-consistent problem within a Hartrcc ap- 
proximation that includes quantum smearing effects from 
the electronic wave functions. Most of these calculations 
are performed in simplified one-dimensional geometries, 
since the self-consistent resolution of the Schrodingcr 
equation becomes prohibitive for an arbitrary disorder 
landscap e) 27 ! 38 . One can hope that our high-field expres- 
sion for the density will turn out to be a very useful tool 
in this context of the study of electron-electron interac- 
tion effects in a disordered system. 

IV. ELECTRON CURRENT DENSITY 



E. Electron-electron interactions and screening 

As a result of electron-electron interactions, the po- 
tential V entering into the previous expressions through 
the Fermi function is possibly very different from the 
bare electrostatic potential (related to confining gates 
and random impurities outside the two-dimensional elec- 
tron gas), and has to be determined self-consistently from 
screening theory. Previous work, rooted in the semiclassi- 
cal picture, used expressions for the electron density such 
as Eq. I|52p as a starting point for Thomas-Fermi type of 



J 



A. General expression 

The local electron current density is defined in terms 
of electronic Green's function by 



j( r > w ) 



eh 
2m* 



(V r / - V r 



m*c 



G<(r, r» 



(67) 

In a first step, this expression can be written in terms of 
vortex Green's functions. Inserting expression (1441) . we 
get 



[£w E [*™,R(r)V r ^ )R (r)-^, iR (r)V r * TOiR (r) + 2i^A^, )R (r)* miR (r) 



+ 00 , ]2 s n. 



k=0 



(68) 



The dependence on the variable r is contained only in the wave functions and the vector potential A, which are all 
known (we remind that we have chosen the symmetrical gauge to write down explicitly the vortex wave functions). 
Using the relation 



's/2l. 



-*m-l 




(69) 



we can rewrite the bracketed term in (l68l) as 



i*7 r* MW Ml v ; _ ^ ( V m ' + 1 *™-R( r )*m'+i.R( r ) - V / i+T$ ra+ i ! R(r)f^. R (r) 
r ^™.RWvl> m , R (r)j - ^ iOTMr)$ . (r) + iVm+T* m+1)K (r)*^ iR (r) 



(70) 
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Inserting expression (|70[) in formula (|68|) , we finally get our starting point for the computation of the current density, 



2m* 



. _ , . n/2 fd 2 R ^ 



'+i.R( r ) _ a/to + l* m+ i jR (r)\E'^, ;R (r) 
»Vm' + l*m,R(r)*^, +liR (r) + Wm + l* m+ i, R (r)*;, R (r) 



+ 00 



71) 



where p(r,u) = — iG < (r, r, w) is the local spectral func- 
tion. Similar to Sec. IIII1 we wish to collect all contribu- 
tions up to order l 2 B to the local current density. 



B. Electronic current at leading order 

The procedure to compute the different contributions 
to the current density is completely analogous to the cal- 
culation of the electronic density done in Sec. IIII1 al- 
though more lengthy. The leading contributions are eas- 
ily seen in Eq. (f7Tj) to come from and 



1. Contribution from : Density-gradient current 

Leading order Green's function is purely diagonal with 
respect to the Landau-level index to, so that we have to 
consider combinations as \Jm + 1 i S rn 1 i!* n+1 in Eq. (fTTj) . 



the strict Is 
sion, 



limit to recover a semiclassical expres- 



Inserting both p^> from Eq. ([50)) and g^ from Eq. 
(considering only the contribution with k = 0) and us- 
ing the useful relation ([54]) . we readily obtain after the 
frequency integral 



j(o) (r) 



-z x V r YW E ^c"-f(6™(R-)) 
El*P.H(r)| 2 -^f^- . (72) 

.p=0 



This contribution to the current density has the property 
that its volume average vanishes: / <i 2 rj(°)(r) = 0. 

As done previously for the local electronic density, the 
density-gradient contribution (|72p can be expanded in 



+00 

j (0) (r)4£ 

rn— 



m 



fitu c V r n F (£ m (r)) x z. (73) 



This result coincides with the formula for the "edge" elec- 
tronic current density derived within a different method 
in the Ref. [36|. It has clearly the form of a current flow 
responding to a gradient of the density. It thus vanishes 
in the incompressible regions where the density is qua- 
siconstant, and becomes important in the compressible 
regions of the system where the local density is strongly 
inhomogeneous and the bare potential is almost perfectly 
screened. Such regions are not necessarily located at the 
edges of the system, but are rather spread throughout 
the system. The denomination of "edge current" is thus 
in some sense abusive. Therefore, we prefer to call con- 
tributions (|?2"j) and ([73")) a density-gradient current. 



2. Contribution from g^ : Drift current 

As emphasized in our previous paper—, the well-known 
drift contribution to the current density appears, in fact, 
beyond the limit lg — > 0, i.e., when considering Green's 
functions g^ 1 ' which take into account the first processes 
of Landau level mixing. Such a drift contribution is how- 
ever of the same order as the density-gradient contribu- 
tion [the reason is that there is a prcfactor Ig 1 in the 
general expression of the current density, see the second 
term in the rhs of Eq. fl7T])] ■ 

Using the general expression for the current density 
(|7T|) with k — and inserting g^ from Eq. (|3"Tj) . we get 



j (1) (r) 



eh . (1) eh fd 2 R n F (£ m+1 (R)) - ra F (£ m (R)) f 

z x V r p<- >(r) + — j-^-^ 



2m 



2ttI 



B m=0 



fkd r 



Im 



+Vm+l\/m + 2( Rc J (d x V + idyV)V m ,ii(r)** m +2,B.( r ) 

I 



(m + 1) |* m+ i. R (r)| 2 z x VrV(R) 

(74) 



In Appendix [C] we provide a detailed calculation of this expression (|74p , which contains a peculiar term (the last 
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one) involving vortex wave functions with Landau indices 
that differ by 2. As a result, we find that the leading 
contribution to Eq. (|74|) reads 

r- +00 

j«(r) = e - d 2 R^2 l*m,R(r)| 2 n F (£„(R))V R y(R)xz. 

m=0 

(75) 

In the limit In — ► 0, above contribution (f75|) yields 

+OG 

j (1) (r) = «F(^(r))V r y(r) x z. (76) 

m=0 

We thus recover the well-known drift current that can 
be found in the literature^, while expression l|75j) consti- 
tutes a quantum version of this drift current, that may 
be used at low temperature. 

C. Electronic current at order 1% 

We aim at collecting exhaustively all contributions to 
the current density that are proportional to l 2 B : This is 
the order where the first dissipative features are expected 
to appear (see Sec. IVI[) . We will simply list here the 
various origins for these terms and refer the reader to 
Appendix [C] for the detailed calculation. 

There is first a subdominant contribution coming from 
g( x \ with expression (|T4[) . In Appendix [Cl this contri- 
bution is denoted j^ b , and is given in Eq. (|C6[) . An- 
other contribution of order l B arises with second-order 
Green's function g^ 2 \ The latter contains diagonal ele- 
ments (to = to') that combine with the first and second 



+00 

j(r) = ^ M£m(r)) 

m=0 



This expression constitutes one of the main physical re- 
sults of the paper and is further analyzed in Sec. [Vl] deal- 
ing with the noncquilibrium transport properties. 



V. CHECKING OUR THEORY: COMPARISONS 
WITH AN EXACTLY SOLVABLE CASE 

A. One-dimensional parabolic confinement model 

The aim of this section is to benchmark our results 
for the local equilibrium charge and current densities ob- 
tained with the vortex states. For this purpose a com- 
parison to the exact solution that can be obtained for 



terms in the rhs of Eq. (fTTj) . and also off-diagonal ele- 
ments 5 mm r±2, which have to be inserted in the second 
term of the rhs of Eq. ([7Tj) . The final expression for 
jw, which also includes the contribution from the func- 
tion g(°^ appearing with the term k = 1 in Eq. (|7ip . 
is given by Eq. (|C8j) . Finally, the off-diagonal elements 
S m ,m'±i of g( 3 ' , calculated in Appendix [Aj combine with 
the contribution from the function g^ associated with 
the term k = 1 in Eq. (fTTj) . giving the final result for j( 3 ) 
in Eq. dCTOl) . 



D. Semiclassical current: The strict Ib — > 
expansion 

As previously done with the local electronic density, 
it is also possible to express the current density under 
the form of a strict expansion in powers of the magnetic 
length. In this section, we want to obtain the correc- 
tions of order l 2 B to the well-known semiclassical expres- 
sion (|76|) for the drift current, which is purely transverse. 
All these sublcading contributions are collected in Ap- 
pendix [C] Since the semiclassical expansion in Ib is only 
valid in a "high" temperature regime, we will only present 
here the terms proportional to the Fermi factor that are 
dominant in this regime with respect to the other terms 
involving derivatives of the Fermi factor. 

Collecting Eq. (|76|) with the contributions from 
Eqs. (|C11[) - (|C14[) , we get the leading contribution to the 
semiclassical current, 



x z. (77) 

I 

the case of a one-dimensional parabolic confining poten- 
tial turns out to be quite cntlightening. We will therefore 
focus here on the following potential profile: 



V(x) = ^m*u; 2 x 2 . (78) 



Following Ref. HH, the exact eigenstates and eigenen- 
crgics corresponding to this particular choice of one- 
dimensional potential are given in the Landau gauge 



V r V + l 



(V r V ■ V, 

fkU r 



■V r V + -li 



AV r V 
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A = Bxy by 



V2"+ 1 n! 7 r 3 / 2 L 



E 



up 



U+ 2 



(79) 
(80) 



obtained from the vortex calculation, consists in devel- 
opping both Eqs. (f8Tj) and (|83| in a strict magnetic length 
expansion at order i^. For this purpose, we develop fl 
and L using Wq/w^ = l\V" (x)/hu> c , and perform a Tay- 
lor expansion. The resulting Gaussian integrals are com- 
puted using the formula 



where O = y 1 u 2 + u>q and L = h/m*Q arc the renor- 
malized cyclotron pulsation and magnetic length respec- 
tively, and H n denotes the n th Hermite polynomial. 
These wave functions, fully extended plane waves along 
constant energy contour while strongly localized in the 
transverse x direction, are certainly very different from 
the vortex states [Eq. ©J, which arc localized in all direc- 
tions without preferred symmetry, so that this compari- 
son provides a very stringent test on the vortex theory. 

Physical observables, such as the local electronic den- 
sity, are readily obtained as 

+°° f+OO 

p(x) = J2 dp |*„ P (r)| 2 n F (E np ) , (81) 

n=0 J -°° 

while the equilibrium current density, directed in the y 
direction, reads^: 



j(x) = ^2 dpj np (x)n F (E np ), 



(82) 



where 



\e\Q 2 ruj 



uj c v n 

+ \e\u c l 



-£pL*-x) \Vnp(r)\' 
*n P (r)\ 2 . 



2 V'ix) 
3 TlhJ r 



(83) 



B. Checking analytically the semiclassical 
expansion 



J —oo 



m n-\ — 



(84) 



This leads to the result 



p(x) 



+oo 
B n=0 



n F (x) + l B n F (x) 



V"(x) 
hio c 

(v'(x)f 



+l%n' F (x)V"(x) + l fn' F (x) 
+| (n + l)n F ( X )(V'(x)f 



(85) 



where n F (x) = n F [(n + \/2)hu) c + V(x)], which is obvi- 
ously equivalent to Eq. (|64l) . 



The calculation of the semiclassical current density at 
order l 2 B follows the same lines, and using formula 



/ d£e"« £ 4 if 2 (0 = 3^2' 1 - 1 n! (n 2 + n + 1/2) , 

J — OO 

(86) 

we recover the leading density-gradient and drift contri- 
butions, 



i i +°° 

n=0 



n F {x)V 



n + - ) hhj c ri F {x)V' 



(87) 



The first check, which is crucial for demonstrating the 
mathematical consistency of our semiclassical limit, as 



while the terms of order l 2 B read 



i, 



ll £ ( n2 + 71 + V2) ^c(V'Y + n' F {x) (\n 2 + \n + ^\ V'V"fiw c + n' F (x) 



2huo, 

+■ I) (VT + n' F (x) 7 - L + 1) V'V" + n F {x)V'^ 



l\3 



r 



One can easily check that collecting all terms in expres- method to generate the semiclassical expansion. We em- 
sions (|C11|) - (|C14I) for a one-dimensional potential yields phasize however that our semiclassical results were de- 
the same result, giving a good confidence in the vortex rived for an arbitrary two-dimensional potential, so that 
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extra terms are actually present in our semiclassical for- 
mula with respect to Eq. ([55]) . In particular some of the 
terms appearing in Eq. (|77p involves derivatives of the 
potential in two orthogonal directions and cannot be in- 
fered from this simple calculation of a one-dimensional 
parabolic potential. It is interesting to note that in 
noncquilibrium it is precisely these additional terms in 
the current density that are seemingly associated with 
dissipative features (see Sec. IVI|) . 



C. Systematic numerical comparison to the vortex 
theory 

We aim here at giving a more quantitative comparison 
for the electronic density between the exactly solvable 
model and the various expansions discussed in Sec. IIIII 
We therefore compute numerically the expression l|8ip 
and investigate both the semiclassical approximation (15^)) 
and the order l 2 B quantum expressions (|50|) and (|58|) . 
Note that the l B corrections [Eq. (|64|) ] to the semiclassical 
result present derivatives of the Fermi factor, which are 
either very small (at high temperatures) or very singular 
at low temperatures, so that they are not included in the 

to 



comparison. In contrast, the l B corrections [Eq. 
the leading quantum result (|50|) are more regular, and 
their inclusion is important to reach quantitative agree- 
ment in an intermediate temperature range, as we will 
demonstrate below. 

In order to show that the reliability of these different 
approximation schemes are rooted in specific tempera- 
ture regimes, we present results for different tempera- 
tures at a given confinement energy hojQ = w c /5, small 
enough to ensure the smoothness of the external poten- 
tial, yet already sufficiently large so that the semiclassical 
approximation is in trouble at low temperature. The elec- 
trochemical potential is also fixed by taking fi* = 3oj c , so 
that three Landau levels are present at the center of the 
system. 

At temperatures not too low compared to the cyclotron 
frequency (first panel (a) of Fig. Q] for fc^T = Q.2hui c ), 
the semiclassical result is still close to the exact solu- 
tion, exactly matched by the quantum result. Lower- 
ing further the temperature (second panel (b) of Fig. Q] 
for fcsT = 0.1hiu c ) shows increasing deviations with the 
semiclassical result, while the complicated variations in 
the exact density are perfectly reproduced by the quan- 
tum formula. In particular, both the small compressibil- 
ity in the filling factor n = 3 plateau [the density of the 
third filled Landau level is slightly greater than the value 
3 x (2ttI b )^ 1 ] and the broad smearing of the n = 1 and 
n = 2 plateaus are quantitatively described. In the very 
low temperature regime, small shoulders appear at frac- 
tional densities (panel (c) of Fig. Q]for fc^T = 0.01Sw c ), 
which are associated with the zeros of the Hermite poly- 
nomials in Eq. (|79|) . These variations are only partially 
reproduced by the quantum expression, but the overall 
agreement remains very good. 
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FIG. 1: (Color online) Local electronic density p(x) in units 
of 1/2ttZ| for the one- dimensional parabolic potential with 
ujo = cj c /5 and n* = 3u> c , as a function of x/Ib, comparing 
exact expression (|8ip (solid curve, label Ex.) with the semi- 
classical expansion (|52|l (dotted-dashed curved, label Sc.) and 
the quantum expansion given by Eqs. (|50[) . (|56|l and (|58[) 
(dashed curve, label Qu.). The three different panels (a)-(c) 
correspond to decreasing temperatures ksT jhwc = 0.2, 0.1, 
and 0.01. 
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D. Zero temperature limit: Resummation of the 
quantum development 



Wc finally motivate the need for a resummation of the 
quantum expression to arbitrary order in Ib in the very 
low temperature regime, as hinted in Sec. Ill Ft As the 
scmiclassical expression for the electronic density (|65l) is 
clearly divergent at low temperature, one can indeed ask 
whether the leading quantum result (|50[) and its order l B 
corrections [Eqs. (|56[) and ([58)) ] give satisfactory results 
for all temperatures. Despite the excellent agreement ob- 
served above, the Fermi factor derivatives appearing in 
these order l 2 B terms tend to give important and uncon- 
trolled contributions in the zero-temperature limit. To 
see this, let us forget for the time being the (negligible) 
terms inversely proportional to hw c in the quantum ex- 
pression for the density, which then simply reads 



p(r) 



r £|* m ,R(r)| 2 ^&n(R)) 



B m = 



m=0 



2ttZ 

I 2 

^-A R n F (£ m (R)) 



p(°)(r)-^A r p(°)(r). 



(89) 
(90) 



Here we have made an integration by parts to rewrite 
the second term in the rhs of Eq. (f8"9"|) . Because p' ' can- 
not change on a scale smaller than Ib, as is clear from 
Eq. ([50]) , the l 2 B correcting term in the rhs of Eq. (|9"0|) 
cannot become singular in the zero-temperature limit, 
in contrast to the semiclassical expression ([6"5"|) . How- 
ever, p' ' does change on the scale Is at the boundary 
of an incompressible region at very low temperature, so 
that the correction becomes of order one and needs to 
be resummed to all orders. The need for a resummation 
is mathematically related to the fact that nonlocal vor- 
tex Green's function has been developped at coinciding 
points in Eq. (|22p , while keeping a finite number of contri- 
butions. A clear example of such a nonlocal resummation 
to all orders is the relation (|44|) between vortex and elec- 
tron propagators. In fact, the correction in Eq. (jiTO)) is 
the combination of the — (^/2)Ar(/ ) term in Eq. (j4"4")l 
and the (^/4)Ar^°) contribution that can be extracted 
from in Eq. ([3%)) . 

By inspecting the recursion relation (|26p in the small 
Ib limit, it is possible to infer that this class of most 
singular terms in the vortex propagator is given to all 
orders by 



1 fl 



k=0 



^(r) = W'Eh( t Ar ) 5 ™ 0) ™ (r) (91) 



so that their combination with Eq. (|4"4")) leads to the fi- 
nal quantum expression for the density in the small but 



nonzero Ib limit 



.R 



fc=0 



^-l-^A R j n F (Um- (92) 



Using an integration by parts, this equation can be writ- 
ten in the equivalent form 



>(r) 



d 2 R , 

m—0 



^ M 

k=0 



n F (£ m (R)) 



4 A k) l*m,a(r)r. (93) 



We note from Fourier analysis that the differential oper- 
ator 



k=0 



^A 



R 



is nothing else than a convolution operator with the ker- 
nel e~ u2 / 4t /(47rf), where t = We can apply this 
to the vortex density, to find 



i-mto +°° 

p B ^°°(r) = /— T ^n F (6„(R))|$ m ,R(r)| 2 (94) 

l$ m ,R(r)| 2 = |^ e - u2 / 4t |* m ,R-u(r)| 2 . (95) 

Performing formally the remaining Gaussian integral 
over u in Eq. (|95p . we find 



(r) 



i g np (um Am(R 



2nl% 

B m =0 



irmU B 



x exp 



(R-r) 
1% 



where A m is the following polynomial: 



Am(R) 



ds r ' 



1 



1 



■ exp 



R 2 2s 
l%TT t 



(96) 



(97) 



s=0 



Final expression (|9"6"|) for the density is applicable down to 
zero temperature, and provides the leading contribution 
in the small Is limit. In the case of a one-dimensional 
potential V(x), it is easy to check that the Gaussian in- 
tegral over the coordinate Y in Eq. (|96[) leads to the 
expected Hermite polynomials. Regarding the remaining 
contributions that can be gathered from Eqs. (|56| and 
(|58|) . involving Landau level mixing processes, a com- 
plete resummation scheme amounts to extra shifts in the 
energies, as discussed in Sec. IIIIDI A final comparison 
is given in Fig. [21 which shows that, as long as loq <C lu c , 
these improved quantum expressions arc undistinguish- 
ablc from the exact result. 
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FIG. 2: (Color online) Local electronic density p(x) for the 
same parameters as in Fig. [He) , focusing on the compressible 
region between filling factors n — 2 and n = 3. This compares 
exact expression (|8ip (thick solid curve, label Ex.), undistin- 
guishable from quantum expression (96) resummed to infinite 
order in Ib (thin solid curve, label Qu"), with semiclassi- 
cal expansion (|52|l (dotted-dashed curve, label Sc.) and the 
quantum expansion up to order 1% given by Eqs. (|50|l . (|56[) 
and (|58p (dashed curve, label Qu.). 

VI. NONEQUILIBRIUM PROPERTIES 
A. Distribution function and irreversibility 

We have solved so far Hamiltonian (O within the 
high magnetic field expansion without fully specifying 
the potential-energy term V(r). This scheme allows us 
to study the equilibrium and noncquilibrium situations 
on an equal footing. At equilibrium, V(r) consists of 
a fixed background potential (including a confinement 
potential and an impurity random potential) and of a 
Hartree potential resulting from the mutual Coulomb in- 
teractions between the electrons. As a result of a self- 
consistent calculation, this yields a global effective elec- 
trostatic potential V e s associated with local microscopic 
electric fields. In the noncquilibrium case, there is in ad- 
dition an external potential-energy contribution reflect- 
ing the appearance of macroscopic electric fields and of 
macroscopic chemical-potential gradients in the system 
induced by the presence of a macroscopic current flow. 
Within the noncquilibrium regime, which is considered 
from now on in this section, the potential term V(r) in 
the Hamiltonian ([5]) consists thus of two different parts, 

V(r) = ^ cff (r)+e$(r). (98) 

Here is the nonequilibrium electrochemical potential 
that now varies in space. The latter term takes into ac- 
count the presence of a macroscopic electromotive field 
E = -V r $. 

In Rcf. [35| we have solved the equation of motion 
for correlation Green's function G < in the vortex rep- 
resentation using the high magnetic field expansion, and 



have established that latter Green's function expressed in 
the vortex variables is related to retarded and advanced 
Green's functions at any order of our expansion in the 
nonequilibrium stationary regime as 

- iG< ;i/ » = in F (fj) - G£ ;B ») . (99) 

In the electronic representation, the quantity —iG < 
which has the character of a distribution function thus 
becomes 

-»G<(r,r» = ^*: 2 (r')vMr)m F (w) 

x(G*.,>) (100) 

where we have used Eq. (f9"9"| . The fact that the same 
relation [Eq. (|100[) ] holds in high magnetic fields in 
the equilibrium regime as well as in the nonequilibrium 
stationary regime can be understood as the realization 
of a local hydrodynamic equilibrium (or quasiequilib- 
rium). This result, which has been established from the 
microscopic derivation of the quantum kinetic equation 
(see Ref. |H), is physically expected given that the mi- 
croscopic characteristic lengthscale for the electron gas, 
namely, Ib oc i?" 1 / 2 , becomes in high magnetic fields 
the shortest lengthscale. This means that it is possi- 
ble to divide the system within a continuum description 
into elementary subsystems which are almost isolated 
from each other, permitting the introduction of thermo- 
dynamic variables depending on the space variable r (see, 
e.g., Ref. [H and references therein). 

It is worth mentioning that we have not introduced so 
far in the resolution of the Dyson equation any averag- 
ing to account for the presence of a random potential. 
In fact, this is not needed at this level since within the 
high magnetic field expansion all the physics in the vortex 
representation appears to be purely local, the Hamilto- 
nian being diagonalized in a closed-form order by order 
in powers of the magnetic length with the use of local 
vortex Green's functions <7„ li:m2 (R). Dissipation and ir- 
reversibility, that are usually introduced already at the 
level of the Dyson equation with the impurity averaging 
procedure to account for the presence of random scat- 
tering interactions in zero or weak magnetic fields, take 
its roots within a different mechanism in high magnetic 
fields. 

In fact, the stochastic character is intrinsic to our high 
magnetic field expansion making use of the vortex basis, 
and we can associate somehow the transformation from 
the vortex to the electronic representations with a loss 
of information (thus irreversibility) provided that there 
exists some dynamical instability in the system. Indeed, 
by solving the Dyson equation in the vortex represen- 
tation, we have basically augmented the set of allowed 
quantum numbers since the vortex basis is overcomplete. 
Note however that the expansion of the matrix elements 
of the potential in the vortex representation is granted in 
high fields due its unicity, which clearly results from the 
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possibility to truncate the series expansion in the mag- 
netic length I b , see Ref . Coming back to the electron 
representation, an indeterminacy illustrated by the pres- 
ence of weight factors, namely, the wave functions ^(r) 
in Eq. (|100p . appears, giving rise to a statistical-like de- 
scription. Our high magnetic field expansion alone con- 
tains thus microscopically the stochastic character which 
is a prerequisite ingredient for the expression of a loss 
of information. This loss becomes effective as soon as 
a dynamical instability associated with a divergence of 
neighboring trajectories exists in the system. This gen- 
eral view is corroborated by explicit calculations in the 
following. In Sec. IVI CI we shall derive a microscopic 
expression for the conductivity tensor that is indeed as- 
sociated in the noncquilibrium regime with an electro- 
chemical potential drop occuring only in the vicinity of 
a local instability of the dynamics which, in the present 
case, is brought by the presence of saddle-points of the 
local equilibrium density. 

It is worth noting finally that we do not take into ac- 
count the interaction of the system with an external en- 
vironment. This alternative approach to dissipation con- 
sidering the system plus reservoir couples the relevant 
quantum system to a large number of environmental de- 
grees of freedom, such as phonons or quantum fluctua- 
tions of the electromagnetic field. Dissipation arises then 
because the system of interest can exchange energy with 
the rest of the larger system. More precisely, a loss of 
information is usually explicitly introduced in the calcu- 
lations when tracing out these environmental degrees of 
freedom. We do not consider this mechanism of dissipa- 
tion as being the most relevant here. On the contrary, the 
irreversibility mechanism described in this paper takes 
place in the bulk of the system iself, i.e. in the two- 
dimensional electron gas. This is a fundamental and key 
point in our transport theory. The relevance of a given 
approach to dissipation can finally be appreciated at the 
level of the comparison between theory and experiments, 
since the dissipative transport properties are strongly and 
intrinsically related to its dissipation mechanisms. 



B. Nonequilibrium current density 

Since the two potential-energy terms V e s and e<I> in Eq. 
can be treated technically on an equal footing in our 
high magnetic field theory, the noncquilibrium current 8 j 
can, in fact, be rather straightforwardly deduced from the 
expressions of the equilibrium current density derived in 
the former Sec. IIVI and Appendix [Cl 

The primary goal of this paper is not to provide a full 
quantitative analysis of macroscopic transport properties 
but just to show that our theory does contain information 
on microscopic dissipative mechanisms, and thus allows 
us to fully determine the spatial dependence of the elec- 
trochemical potential <f>(r). For the sake of simplicity, we 
shall therefore restrict ourselves to the regime where the 
current density can be expressed in a local form [local ex- 



pressions ([73]). (|76|) and (|C11|) - (|C14[) ] . This regime does 
not correspond to the lowest temperatures for which the 
nonlocal nature of the current density associated with 
quantum tunneling becomes predominant . 

It is clear, e.g., from formulas (f73|) and (f76|) . that a 
nonequilibrium current can be generated in the linear re- 
sponse and at a uniform temperature by simultaneous 
density and electrostatic potential variations. This indi- 
cates that in principle we can separate the total current 
into two different contributions. One contribution corre- 
sponds to the diffusion current (terms involving gradients 
of the density) whose physical origin is associated with 
the tendency of the system to make the density uniform. 
The other contribution represents the current produced 
by the electric field which accelerates the electrons: It 
corresponds to electrical conduction which occurs by def- 
inition for a uniform density. The true driving force for 
the electrons is finally a combination of chemical and 
electrostatic potentials differences, i.e., it is character- 
ized by the variations of the electrochemical potential 
4> which yield at the macroscopic scale a voltage drop. 
Although the roles of the electrostatic and chemical po- 
tentials are fundamentally different microscopically, the 
precise composition of $ is irrelevant in the linear re- 
sponse regime. Indeed, the integrated nonequilibrium 
electrical current can be seen, e.g., either as resulting 
entirely from a density-gradient current, or equivalcntly, 
as being entirely produced by macroscopic electrostatic 
variations (note that this equivalence is only valid in the 
linear response). In this paper, we shall adopt the point 
of view of the conduction mechanism, i.e., the case where 
a transport current <5j is only sustained by electrostatic 
variations (the electrochemical potential changes arc only 
identified with the electrostatic potential changes in the 
system). 

Any analysis (semiclassical or quantum) of the 
noncquilibrium properties where interaction effects are 
expected to play a crucial role involves the simultane- 
ous resolution of a transport equation and of the Poisson 
equation (in this direction, see, e.g., Ref. H^). Since the 
equilibrium and nonequilibrum regimes can be described 
with almost the same expressions, we can first consider 
at a qualitative level that the interaction effects in the 
noncquilibrium case do not differ substantially from that 
known in the equilibrium case. The screening in the com- 
pressible regions being almost perfect at low tempera- 
ture, we expect a priori that the noncquilibrium conduc- 
tion current is principally confined to the incompressible 
regions where most of the macroscopic electrostatic vari- 
ations giving rise to voltage drops can occur (in other 
terms, this means that only one type of currents - diffu- 
sion or conduction - contributes to a given region in the 
ideal case of perfect compressibility and incompressibil- 
ity). This aspect concerning the noncquilibrium conduc- 
tion current distribution has already been put forward 
by different author a 29 i 30 i 40 . 

Using Eqs. O, (J7SJ), OH and (ICTB - (|CT1 . and keep- 
ing only the terms that are linear in variations of the 
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electrochemical potential $ (since we consider the linear 
response) and that do not contain derivatives of the Fermi 
function factor (we consider the conduction mechanism 
which involves the whole Fermi sea), we get at leading 
order for the nonequilibrium conduction current 

e 2 

<5jo(r) = — n F (U(r)) V r $ x z. (101) 

m=0 

I 



2 +°° 

5h(v)= e -Y. n ^{v))l 2 E 

m=0 



From now on the Fermi factor 7ip(£ m (r)) is a functional 
of the effective equilibrium potential V^ff which differs 
only slightly from the bare potential in the incompressible 
regions of the system where the screening is ineffective. 
A smooth spatial variation of the factor np(^ m (r)) exists 
in these regions as a result of the finite temperature. 

Obviously, the leading contribution (|10ip yields local 
Ohm's law which takes the form 

<5j (r) = a(r) E(r) = a H (v) z x E(r), (103) 

with a local conductivity tensor containing only the 
transverse Hall component 

2 +°° 

^(r) = -r £ n F (£ m (r)) , (104) 
m=0 

where £ m (r) = E m + V e s(r). Ohm's law (|103[) with local 
Hall coefficient (|104[) is already well-known and is used 
in most of the existing transport theories discussing the 
integer quantum Hall effect. The absence of diagonal 
components for the conductivity tensor in this (semiclas- 
sical) limit Ib — * is rather welcome, since it is com- 
patible with an extremely small longitudinal resistance 
as observed when the Hall resistance presents plateaus. 
However, this absence points out at the same time an 
insufficiency of the formula (| 103[) to describe the transi- 
tion region between the Hall plateaus when high peaks of 
the longitudinal (dissipative) magnetoresistance are seen. 
This insufficiency is, in fact, cured when considering con- 
tribution (|102[) to the current arising from the next order 
terms in the Ib expansion, as shown further. At a general 
level, we note that our quantum-mechanical derivation of 
the transport current justifies on a microscopic basis the 
use of phcnomcnological models assuming a local conduc- 
tivity tenso r 28 i 29 i 30 i 31 i 32 i 33 i 41 ' 42 that have been considered 
so far to explain successfully some transport features of 
the quantum Hall effect. 

To our knowledge the first quantum corrections [Eq. 
(I102[) ] to Ohm's law (|103[) had not been derived before 



From Eq. (|77|) . we get a correcting contribution to the 
nonequilibrium current <5j which is second order in Z^, 



x z. (102) 

I 

in the literature. We find that they contain local cor- 
rections (which give rise to both transverse and diagonal 
components in the local conductivity tensor) as well as 
nonlocal corrections (terms involving second- and third- 
order derivatives of $). These nonlocal terms can be 
viewed as fingerprints of the nonlocal quantum tunneling 
processes in the considered semiclassical regime. 

C. Spatial dependence of the electrochemical 
potential 

The expansion of the current density in powers of Ib 
has led us quite naturally to a local continuum descrip- 
tion of current conduction. Within this "classical" pic- 
ture of transport (our theory is nevertheless developed in 
a fully quantum mechanical framework), the stationary 
equation of continuity ensuring the charge conservation 

V r -j = 0, (105) 

supplemented by boundary conditions, constrains the 
spatial dependence of the electrochemical potential when 
applied to the nonequilibrium current density 5] pro- 
vided that some dissipation mechanisms are accounted 
for within the considered expressions for <5j [note that 
Eq. (|105p becomes an identity for the equilibrium cur- 
rent density, as can be easily checked] . 

Inserting the leading contribution 8 jo [Eq. (|103[) ] into 
Eq. (|105p and using E = — V r $, we get the equation 

(V r cr H x V r $) • z = 0, (106) 

which has been thoroughly discussed in the 
literatur o 28 i 31 ' 32 . From this equation, it turns out 
that the electrochemical potential lines and the lines of 
constant an must coincide. It is worth noting that the 
condition (| 1 06[) is automatically obeyed at the critical 
points of <7ff which correspond also to the critical points 
of the density, or of the potential V e g. This means 
that there still exists a degeneracy in the vicinity of 
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these critical points, which has to be lifted. Although 
being small, the correcting contributions [Eq. (|102p ] will 
play this important role of dictating locally the spatial 
dependence of <&(]"), as we prove now. 

Since we are considering the nonequilibrium current 
density in the neighborhood of V r T4ff = 0, we can first 
safely ignore in Eq. (|102p the term proportional to 'V r V c g 
and which involves the second-order derivative of <f>. At a 
preliminary stage, we shall also disregard the other non- 
local term (with the third-order derivative of $), and 
justify this assumption a posteriori. Consequently, the 
second-order contribution to the current reduces to 



Sh(r) 4l M£m(r))Z| (Vr f " Vr) V r y eff x z. 



Combining this Eq. ([TUT]) with Eq. ([TO?]) , we get local 
Ohm's law with a local conductivity tensor being given 
by 

= <Mr) + hf 1 [Vcs] (r) ) ' (108) 

where 1 is the 2x2 identity matrix and H is the Hessian 
matrix of the function V e g, i.e., we have H. [Veg]- (r) = 

We remark that the local conductivity tensor does not 
exhibit the usual symmetries, i.e., the Onsager-Casimir 

I 



where the notation Tr means the trace. In the neighbor- 
hood of a critical point, which for practical convenience 
is taken at the origin (x = y = 0), we have 

V r <r H {r) w (r- V r ) V r a H \ . (110) 

The Hessian matrices of the function o\ff(r) and of the 
function V c g{r) being proportional at the critical point, 
we can choose, without loss of generality according to 
the form of the Eq. p09p . the x and y axes such that 
both Hessian matrices are diagonal. This means, e.g., 
that V e g (r) is expanded close to the origin as 

2 2 

Kff(r) = Veff(0) + ay + 6^- (111) 

where a = d^V g\„ and b = SLKjffL. For this situa- 
tion, Eq. (|109p becomes then 

(a - b) 8L* + 4 [ax8y* - byd x <S>] = 0, (112) 



reciprocity relations. For example, we find here that gen- 
erally a xx {r) = -<J yy {r) and a xy (r) ^ -a yx (r), whereas 
the Onsager relations imply a xx = a yy and a xy = —o yx . 
In fact, it is worth noting that the Onsager relations re- 
sult from fingerprints of the time-reversal invariance of 
the microscopic equations after some averaging proce- 
dure (see, e.g., Ref. l39t ). In the present case, we have 
derived a local conductivity tensor from the microscopic 
equations without resorting to any averaging procedure. 
Obviously, the local terms involving the Hessian matrix 
contribution vanish in the volume average; the Onsager 
relations are then restored. This indicates that the Hes- 
sian matrix terms can be interpreted as a result of lo- 
cal fluctuations only. Let us also emphasize that current 
()107p purely stems from Landau level mixing processes. 
The found sign difference between the two local diago- 
nal components which appears as rather unexpected and 
unconventional could be seen as a reminiscence of the an- 
tisymmetry imposed by the Lorentz force, antisymmetry 
which is usually only exhibited by the Hall components 
(see the conductivity tensor at leading order). Anyway, 
we shall show in the following that the precise form we 
have found for the local conductivity tensor leads to re- 
liable physical results. 



Inserting second-order contribution (|107|) into conti- 
nuity Eq. (|105[) . condition (| 1 06[) is now replaced by the 
differential equation 



(109) 

I 

with 

+oo 

X2 = 1 'bT^ ■ ( 113 ) 

Y^f^c [-n' p (UO))] 

m=0 

To get rid of coefficients a and b in the differential equa- 
tion [Eq. (|112p ], it is useful to introduce the change in 
variables 

x' = ax - (3y (114) 
y' = ax + (3y. (115) 

If the critical point corresponds to a local extremum 
(situation with ab > 0), we can take 




(116) 



(V r a H x V r $) ■ 



-^CTffTr ■ 



\ 1 )w[^ ff ]w[$]} = o, 
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and Eq. (| 1 1 2[) then reduces to 

a y 2 v $ - d 2 x , x ,<S> + ± [x'd y ,$ - y'd x ,$\ = 0, (117) 

where e = +1 if \a\ > \b\ and e = -1 if \a\ < \b\. The 
general solution of Eq. (|117|) is 



$(x',y' 



A + Be 



[C + D(x' 2 +y' 2 )],(U8) 



where the coefficients A, B, C et D are constants of 
integration. As boundary conditions, we require that 
the electrochemical potential tends to constant values 
far from the critical point. We therefore necessarily get 
D = B = 0. There is consequently no macroscopic volt- 
age drop associated with the crossing of a local extremum 
(the special case a = b ^= 0, which can be readily obtained 
from Eq. (| 1 12[) leads to the same result). 

Now, if the critical point corresponds to a saddle-point 
(situation with ab < 0), we can choose 



Q 



-b 1 



= 



b — a' 



(119) 



so that Eq. (|112|) becomes 



- d 2 x , x ,<s> + ^ [y'dy^ - x'd x ,$] = o. (120) 

Looking for a solution with separable spatial depen- 
dences, we find that the only solution of Eq. (|120p rep- 
resents the product of two error step functions in the x' 
and y' directions, 



$(x',y') = 
where 



A + B erf 



erf(x) 



y/2X 



C + Dcii 



dt. 



y 



V2A 



(121) 



(122) 



We observe with the solution (|121j) that far from the 
saddle-point the electrochemical potential tends to dif- 
ferent constant values depending on sectors. This solu- 
tion is consistent with the picture of four different regions 
characterized by four different electrochemical potential 
values with an electrochemical drop resulting from the 
saddle-point crossing, which is a major ingredient in the 
network models that have been developed to describe the 
peaks of the longitudinal conductance in the transition 
reg ime between quantized Hall plateaus (see, e.g., Ref. 
kll ). Thus, our conductivity tensor confirms at a micro- 
scopic level the special role played by the saddle-points 
of the density in the dissipative feature o 17 ! 28 ' 37 ' 41 ' 43 ! 44 . 

Finally, we turn back to the condition of validity of Eq. 
(|109|) which has been established under the assumption 
that the nonlocal term in Eq. p02[) involving the third- 
order derivatives of $ play a negligible role. Clearly, this 



is justified provided that <£> is smooth enough. Consider- 
ing expression (|113p giving the characteristic lengthscale 
A for the spatial variations of the electrochemical poten- 
tial, we note that this assumption appears fully justified 
as long as the function n F (£ m (0)) remains quite small, 
i.e., as long as the saddle-point filling factor is close to 
an integer. Conversely, we can conclude that the nonlo- 
cal term which is associated with quantum tunneling be- 
comes non-negligible at low temperatures when the local 
chemical potential /i(0) = fi* — V c s(0) approaches a Lan- 
dau level. This regime which occurs for a narrow range in 
magnetic fields will be investigated in detail elsewehere. 



VII. CONCLUSION AND PERSPECTIVES 
A. Summary 

In summary, we have developed a systematic high mag- 
netic field expansion, which permits to find in a recursive 
way, order by order in powers of the magnetic length I b , 
Green's functions for the quantum problem of an electron 
confined to a plane and subjected to a slowly-varying 
potential in high magnetic fields. Using this theory, we 
have derived functional quantum expressions for the lo- 
cal equilibrium density distribution and current density 
at the first two leading orders. These expressions which 
contain Landau level mixing processes in a controlled way 
and quantum smearing effects associated with the finite 
extent of the wave function at finite magnetic fields form 
the starting point for future quantitative investigations of 
screening effects at low temperatures in two-dimensional 
disordered Hall liquids. We have checked the accuracy 
of our general functionals against the exact solution of 
a one-dimensional parabolic confining potential, demon- 
strating the controlled character of the theory to get equi- 
librium properties. Furthermore, we have shown that our 
technique gives a natural and systematic access to semi- 
classical expansions in powers of the magnetic length of 
the physical observables. For example, we have been able 
to derive for the first time the semiclassical corrections 
of order l B for the local charge and current densities. 

Moreover, we have proved microscopically that in high 
magnetic fields the electronic system can be described 
within a local hydrodynamic regime and that the elec- 
trical conduction transport takes a quasilocal form. As 
an important result, we have put forward that our ap- 
proximation scheme with the Is expansion intrinsically 
captures dissipation mechanisms at the microscopic level 
and accounts for quantum tunneling processes. For ex- 
ample, we have derived microscopic expressions for the 
local conductivity tensor, which contains both Hall and 
longitudinal components, the dissipative features appear- 
ing at the order l B , i.e., at finite magnetic fields. Fur- 
thermore, we have established from the special form of 
this local conductivity tensor that a nonzero gradient of 
the electrochemical potential is exclusively generated by 
the saddle-points of the density distribution. A general 
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understanding of the transport properties at the micro- 
scopic level now seems accessible. However, the proce- 
dure of computation of the macroscopic transport coeffi- 
cients, and in particular the treatment of nonlocal effects 
induced by quantum tunneling, requires additional work, 
that is currently under way. 

B. General perspective 

Finally, we want to address a more general perspec- 
tive, which is beyond the quantum Hall effect, namely 
the issue of dissipation in physics and especially in quan- 
tum mechanics. Indeed, we believe that our systematic 
expansion in ascending powers of the magnetic length 
could shed new light on this important issue by illustrat- 
ing how the irreversible evolution of a quantum system 
can emerge from the consideration of microscopic equa- 
tions which are time reversal invariant. The system con- 
sidered in this paper is maybe the simplest system one 
can consider to answer the latter question, since it in- 
volves only 2 degrees of freedom. Interestingly, the mag- 
netic field (which plays the role of a tuning parameter 
here) controls the degree of mixing between these two 
degrees of freedom, which corresponds at the classical 
level to the cyclotron motion and the guiding-center mo- 
tion. In high magnetic fields and for a smooth arbitrary 
potential, this mixing becomes very weak as a result of 
the strongly different timescales associated with the two 
kinds of motion. The system even becomes dynamically 
integrable in the strict limit of infinite magnetic fields 
when the mixing between the two degrees of freedom is 
no more possible. Therefore, the high but finite magnetic 
field regime can be associated with a classical regime of 
soft chaos. At the quantum-mechanical level, it is clear 
that the quantization of the kinetic orbital motion which 
introduces robustness (in the sense that it considerably 
constrains the possible variations of the orbital motion) 
renders this exchange between the 2 degrees of freedom 
even much more ineffective. We thus expect that the 
quantum system is somehow even closer to integrability 



than the classical one. 

In classical chaotic systems (this is, for example, the 
case for the present disordered system in low magnetic 
fields), irreversibility and dissipation are often associated 
with the technical impossibility to fully describe the tra- 
jectories as a result of complicated mixing mechanisms 
between the degrees of freedom. This complexity is then 
transposed in terms of a stochastic description, thus ex- 
pressing a loss of information. We have shown that in 
high magnetic fields it is not required to average over the 
disorder configuration in order to find an analytical ap- 
proximate solution to the quantum problem, contrary to 
the situation at low magnetic fields. Moreover, we have 
noticed that time irreversibility has nevertheless been in- 
troduced at some stage of the derivation, since our high 
magnetic field theory accounts for dissipation features 
related to time-decaying states. It turns out from first 
considerations that the dissipation involves Landau level 
mixing processes and arises from the conjunction of lo- 
cal quantum fluctuations with a local dynamical instabil- 
ity taking place at the saddle-points of the local equilib- 
rium density (a saddle-point is necessarily characterized 
by stable and unstable directions which can be defined in 
an obvious manner). Interestingly, in low magnetic fields 
the electrical conduction is also directly related to an- 
other instability mechanism which is realized by the sen- 
sitivity to the initial condition characterizing the chaotic 
systems. In brief, insights in this general perspective of 
understanding the emergence of dissipation in quantum- 
mechanical systems could be gained from closer investi- 
gations of the high magnetic field expansion developed in 
the present work. 
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APPENDIX A: OFF-DIAGONAL ELEMENTS OF THIRD-ORDER GREEN'S FUNCTIONS 

In this Appendix, we provide the detail for the derivation of the elements i mi , m2 ±i of third-order Green's function 
g( 3 ', that are needed for the computation of the second-order current density performed in Sec. IIVI and Appendix [C] 
The function g^ obeys the equation 



m 3 

+ (d x - idy)V(R)(dx + idy)gM. ma (R) + (d x - i dyK« ;m2 (R)(9x + idy)g^ ima (R). (Al) 
Here is a list of these numerous components S mim ^±i of ffmijmaCR) arising 



from the combination v^g^: 
1 



2u> mi ui m2 



[(mi + l)y/m~[ 5 mi , m2 +i(dx ~ idy) + (m 2 + l)y/m^6 mi +i,m 2 (dx + idy)} AV, 



from the combination v^g^ 
mi + 1 



"mi "'ma 



ArF [v^l <W,m 2 +l(<9x - «<9y) + <5 mi + l,rn 2 (<9x + idy)] V 



+ 



(m 2 + l)y/m2 



2uj 



'ma- l^m 2 w m a +l 



Omi,m 2 + l 

[(dx - ic>y) 2 V] (dx+idy)V 
[(<9x+i9y) 2 y] («9x-zdy)V, 



• from the combination w^ 1 ^*- 2 -': 



m 2 + 1 , m 2 



1 



W m2 + 1 Ulm-2-l 



/m7£mi,m 2 +l(<9x - *<9y) + V^*roi+l,m a (9x + V" 



mi^/?7i 2 



2w mi _iw mi o;„ il+ i 



Om 1 + l,m 2 



2w mi _io; mi w mi +i 
mi^/m 2 

(mi + l)y/m~\ 

• and from the combinations (<9x — idy)v(dx + idy)g: 



- - + 5mijro2+1 + [(dx - idy)V] 



<W, ma +i [(dx + [(dx - idy)vf , 



|v R vf 



1 



1 



J up 1 

a r v 



LU m , OJ 



mi(d x - idy)5 mum2+ i + x /m 2 "(9x + idy)5 mi +i,m 2 ] V 
— {d x - idy)5 mi , m2+ i + ^^-(d x + id Y )6 mi+ i, m2 V 

ii ^m 2 



x/mj 



<W 2+ i [(fibc - idyfV] [(dx + idy)V] 
<W+l,™ 2 [(dx + idyfV] \(d x - idy)V] . 



Regrouping the terms of the same form, different contributions (|A2I) - (|A5[) to the components S mi ,m 2 ±i of 
rearranged as 

• terms with Ar(9x ± idy)V: 



-A R [(mi + l)y/rn^ S mum2+1 (dx - idy) + (m 2 + l)y/m2 (d x + idy)} V/u) mi u> m2 , 

terms with Ar1/(<9x ± idy)V: 
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m 2 + 1 mi + 1 



[y/mi S mu „ l2+1 (d x - idy) + y/m^ S mi +i,m 2 (dx + idy)] V 



/mi 



m 2 



(dx - idy)S mi ,m 2 +i H (dx + idy)8 mi +i,m 2 



V, 



terms with [(d x ± i<9 y ) 2 l/](<9 x T ift-)F: 



(A7) 



m 2 



mi + 1 



^m2-l w ra2 W m 2 +l W mi _i(J mi CJ mi _|_i U> mi U! m2 

mi m 2 + 1 2 



Wmi-l^m^nu + l ^m 2 ~l^m 2 ^m 2 +l UJ mi UJ rn 2 



y/m[S mum2+ i [{d x - id Y ) 2 V] [(d x + id Y )V] 

m2 $mi+l,m 2 

[{dx + idvfV] [(d x -idy)V) 



terms with | V R V\ 2 (d x ± idy)^: 



|v R vf 



m 2 



mi + 1 1 



mi + 1 

' 1 1 2 1 2~ 

WmiWmi + 1 ^m 2 -l^>m 2 U mi LO rn2 U> mi 



m, 2 + 1 



mi 



m 2 + 1 



2 2 

w m2 w m2 +i w mi _iw mi w mi u; m2 w mi w m2 



[(9a- - idy)V] 
m^S mi+lm2 [(d x + id Y )V] 



(A8) 



(A9) 



APPENDIX B: PROOF OF USEFUL RELATIONS 



In this appendix we prove identities (|54|) . and (|Cip - (|C2j) . First, with the help of Eq. (|69|) we can find that 

V2 



(ft - id y ) |* p , R (r)| 2 - j- (Vp*;, R (r)*p_i, R (r) - *;+i >R (r)* P , R (r) 



(Bl) 



what defines in a recursive way the combination s/p + 1 ^t +i R (r)\I/p iR (r). From this relation (|B1[) . it is then straight- 
forward to obtain identity (|54|) . 

Using Eq. (|69|) . it can be readily established that 



^/2 



(ft. - ift) {*«,a(r)*^ +liR (r)} = — [^VwW^i,hW - \/m + 2 * TO)R (r)^ +2jR (r)] 



From this relation (IB2I). we deduce that 



2*«+a,K(r)*m,R(r) = V^*;„+i, R (r)* m _i, R (r) - -^L (ft - «ft) |*;„+i, R (r)* m , R (r)} . 



Multiplying Eq. (|B3[) by V m + 1, we get a recursive relation which yields 

Is 

[o x - ia Vl 

p=0 



\/m + 1 sfm + 



Finally, using identity we get the result Eq. (jCl 
From Eq. (|69|). we can get 



(B2) 



(B3) 



(B4) 



A r |* p , R (r)| 2 = - {(p+ 1) |* p+ i, R (r)| 2 +p|* p _i, R (r)| 2 - (2p+ 1) |* p , R (r)| 2 } . 
Therefore, we can write 



(B5) 



^(m + l-p)A r |f PiR (r)| 5 

p=0 



7Ti 

E( m + 1 ~P) {(P + 1) I*p+i,rW| 2 + P l* P -i.R(r)| 2 - (2p + 1) |* P , R (r)| 2 } 



Z 2 - 

l B p= o 



;2 



(m + l)|* m+ i, R (r)| 2 -£|^ R (r) 



(B6) 



p=0 



which proves identity (IC2[) . 
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APPENDIX C: CALCULATION OF THE ELECTRONIC CURRENT AT ORDER 1% 

In this appendix we present the detailed derivation of the quantum (Appendix IC 1[) and semiclassical (Appendix 
IC 2[) expressions for the electronic current density up to order l^. 

1. Quantum expressions for the current 

a. Contribution from g^ 1 ' 

First-order vortex Green's function has the total contribution to the current given by formula (|74|) . from which 
the leading-order term was extracted in Eq. (|75"|) . We express here the formula ([71]) in a form that makes explicit its 
leading and subdominant contributions. Using the identities proven in Appendix [B] 



Vm + lVm + 2 2 R (r)* TOjR (r) 



7 1 2 m 

IB 



^ ( 5 » - id v) 



£(m + l-p)|vI/ p , R (r)| 2 , (CI) 



p=0 
2 



f A r ^(m+l-p)|vI/ p , R (r)| 2 = (m+l)|* m+liR (r)| 2 -^|vI/ p , R (r)| 2 , (C2) 

p=0 p=0 

we can rewrite the combination ^ m ,Ji( r )^m+2 R ( r ) of vortex wave functions appearing in Eq. (|74[) in the following 
way: 

V^TTV^+2 ^ ^ (d* V + iarV)* m ,R(r)*^ +2>R (r) - l%x x (V R F ■ V r ) V r J2(m + 1 - p) |<J P , R (r)| 2 

+z x V R V (jr |* p , R (r)| 2 - (m + 1) |tf m +i, R (r)| 2 j . (C3) 
\p=o / 

Inserting expressions ([55)1 and (|C3[) in Eq. (|74[) . we then express the current density as 

jW(r) = 5 £ n f (U + i(R)^-n f K m (R)) ^|^ R(r) |^ x VRy(R) 

° m— \p— 

m \ 

+Z|zx [V R V(R)- V r ] V r ^(m + l/2-p)|* p , R (r)| 2 , (C4) 

p=o / 



v r ( V R U • V r |^ P , R (r)| 2 ) = [V R U • V r ] V r ( |1< P , R (r)| 2 ) . (C5) 



where we have used 



After a straightforward simplification of expression (|C4[) . the current density j*- 1 -* can finally be divided into a leading 
contribution given by Eq. (|75[) and a subdominant contribution which reads 



« +oo / m 

j«(r) = t W £ n F (£ m (R)) I| [V R V(R) ■ V r ] V r £ |* p , R (r)| 

•* m=0 \p=0 



2 l*m,R(r)| 2N 



x z. (C6) 



6. Contribution from g^ 2 ' 

Second-order vortex Green's function (|35() contains diagonal elements (m = m'), which contribute within the first 
and second terms in the rhs of Eq. ([71]) . It has also off-diagonal elements <5 mjm '±2, which combine with the second 
term of the rhs of Eq. (|71[) . to give terms involving wave functions with adjacent Landau levels (wave functions with a 
Landau index difference of 3 are also obtained, but these contribute to the current at order lg and will be discarded). 
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After inspection, the contribution from the function </°) appearing with the term k = 1 in Eq. (|7ip combine very 
naturally with these terms from g( 2 > , so that the starting expression reads: 



j (2, (r) 



eh 



1 y ' m* / 9-tt / 9-71-/2 ./o 



"S m=0 



Re 

111! 



{*m+2,a( r )*Tn+i,R(r)flffi+a(^w) 



2m* w m* J2tt J 2ttI{ 

tj m— u * 

+* m+ i,R(r)*;„, R (r) [ 5l (2) „1(R,^) - ARflW<(R,o;)] } . (C7) 
After using Eqs. (|54[) and (|58|) . and performing the remaining energy integration, we can rewrite expression (|C7j) as 



j (2) (r) 



m=0 



f K(£ m (R)) 



mA R V 



[to?i f (£ m _i(R)) 



-(m + l)n F (£ m +i(R)) - (2m + l)n F (£ m (R))] J fi^ c V r 
(m + 1) | [n F (W2(R)) + n F (£ OT (R)) - 2n F (£m+i(R))] 

k(W 2 (R))-» f fc(R))] 



2 |*m,R(r)|' 



El*^(r)l 

.P=0 

v R y(v R y-v r ) |v R vf 



2hw r 



^ m+1 ~| 

E iwoi 2 x *• 

J p=0 J 



(C8) 



c. Contribution from g 



(3) 



Finally, there exist second-order contributions to the current density coming from the elements of third-order 
Green's function g^ which couple adjacent Landau levels [this contribution arises from the second term in the rhs 
of Eq. (|7ip ]. Similar to the previous calculation, these recombinc nicely with the contribution from the function g^ 1 ' 
associated with the term k = 1 in Eq. (JTTJ) . Our starting expression thus reads 



j(3)( r ) = j^L f d2R ^ 



, fv^l^RWI 2 ' 
m* / 2n I 2-kLu '-^ 2 

B m=0 



Re 



(3)< 



(C9) 



Inserting the explicit expressions for the first- and third-order Green's functions [expressions (f5T|) and (|A6[) - (|A9[) ] and 
performing the integration over the energy u>, we finally find after tedious calculations 



j (3) W 



Jd 2 R E ym l*™,R(r)| 2 [n F (£ ro (R)) - n F (e m -i(R))] A R V R V 



(m + 2)n' F (U(&)) -mn' f (^_i(R)) + — [n F (£ ra _i(R)) - n F (£ TO (R))] 

tlUJ r 



a r WrV 



+A{[n F (e ro _i(R)) - n F (£ m (R))] V R V} + [ArVVrF - V R {| V R vf } 
x [4fiw c n' F (f m _ 1 (R)) + (m + 1) [n F (£ ro _i(R)) - n F (U+i(R))] + (m - 1) [n F (£ m (R)) - »jr(^ m _ a (R))]] 
(m + l)7i F (6„+i(R)) - (m - l)n F (£ m _ 2 (R)) + (3m - l)n F (fm-i(R)) - (3m + l)n F (£ m (R)) 



(^ c ) 2 [<(£ m (R)) - n F (^ m _!(R))] + 2/^ c [n F (£ m _x(R)) - n F (£ m (R))] 



2(^ c ) 2 



(CIO) 



2. Semiclassical expressions for the current 



The second-order contributions involve several terms according to their different possible origins. A first term comes 
with the expansion of the density-gradient contribution j' ' [Eq. (|72p ] 



j(0) (r) 



ell ^ (m+lf 



, 2 +oo 
_ l _B \ ~* 

ft 2 ^ 

m=0 



frw c A r V r [n F (£ TO (r))] x z, 



(Cll) 
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where we have used Eq. (|60[) . After making two integration by parts and using Eq. (|60p . the second-order term for 
the current density arising from Green's function [Eqs. (|75|) and (|CJ6[) ] takes the form 



J (1) (r) = ^(3m + 2)A r [n F (^ n (r))V r V] x z. 



/2 +oo 



m=0 



(C12) 



Second-order terms brought by the contribution j( 2 ) are written as 



p I 2 + °° 



m + - ) ?kj c V r 



"F(Cm(r)) 771 A r F 



IlUJr 



\VtV\' 



X z 



-t— K(U(r)) (V r y • V r ) VrV + (V r ■ {n F (Cm (r) ) V r F}) V r F] X Z 



m + - ) [A r W r n F (fm(r)) - n F (£ ro (r))A r V r y - 2 (V r n F (£ m (r)) • V r ) V r V] x z 



(C13) 



Finally, the terms originating from the contribution j' 3 ) yield the following second-order correction to the current 
density 



j (3) (r) = ^|^zxj( m + l)a ra W) 



I 2 + °° 



A r VV r V - V, 



(|v r yf 



n ' Fi t ir)) ~ ) 
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